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ABSTRACT 

A  mathematical  model  developed  for  digital  gunfire  control  using  the 
NSWC/DL  effective  drag  functions  was  modified.  The  NSWC/DL  functions  were 
replaced  with  fourth-order  Runge-Kutta  numerical  integration  of  particle 
trajectory  equations  for  calculation  of  gun  elevation  angle,  projectile 
deflection,  and  projectile  time  of  flight.  Projectile  drag  coefficient 
is  obtained  during  the  integration  process  by  table  look-up  of  prestored 
values.  The  sequence  of  computations  is  given  for  the  present  position, 
prediction,  ballistics  and  gun  order  sections  of  the  model. 

The  current  version  of  the  modified  model  was  designed  not  for  use  in 
a  gunfire  control  system  but  as  a  tool  for  evaluation  of  gunfire  control 
models. 
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SYMBOLS 

FORTRAN 


Model 

Text 

AM 

AM 

Parameter  for  each  projectile  relating  change  of 
Initial  velocity  to  variation  of  weight  from 
standard 

AX.AY.AZ 

AX,Af,AZ 

Components  of  earth's  rotation  vector  in  X,  Y, 
and  Z  directions,  respectively 

B 

B 

Angle  between  ship's  centerline  and  line  of  sight 
measured  in  horizontal  plane  clockwise  from 
ship's  centerline 

BD 

Bd 

Angle  between  ship's  centerline  and  line  of 
sight  measured  in  deck  plane  clockwise  from 
ship's  centerline 

BDGP 

Bdg' 

Angle  between  vertical  plane  through  ship's 
centerline  and  normal  plane  through  line  of  fire 
measured  in  deck  plane 

BG 

Bg 

Angle  between  ship's  centerline  and  line  of  fire 
measured  in  horizontal  plane  clockwise  from 
ship's  centerline 

BGY 

Bgy 

Azimuth  of  line  of  fire  measured  clockwise  from 
north 

BWY 

Bwy 

Azimuth  of  direction  from  which  true  wind  is 
blowing,  measured  clockwise  from  north 

BY 

By 

Azimuth  of  line  of  sight  in  horizontal  plane 
measured  clockwise  from  north 

BYTG 

ByTg 

Azimuth  of  target  at  end  of  dead  time 

BY2 

By2 

Azimuth  of  future  position  measured  in  horizontal 
plane  clockwise  from  north 

CANS 

CANS 

Drift  constant  for  projectile 

CM 

CM 

Mach  number  of  projectile  at  altitude  Z 

CO 

Co 

Azimuth  of  own  ship's  course  measured  clockwise 
from  north 
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SYMBOLS  (Continued) 


FORTRAN 

Model 

Text 

CQO 

Cqo 

Azimuth  of  own  ship's  heading  measured  clockwise 
from  north 

CS 

CS 

Velocity  of  sound  at  altitude  Z 

CSS 

CSS 

Velocity  of  sound  at  standard  atmospheric 
temperature 

DCQO 

DCqo 

Rate  of  change  of  ship's  heading 

DDS 

DDS 

Correction  of  ballistic  air  density  for  non¬ 
standard  projectile  weight 

DEIO 

DEio 

Rate  of  change  of  pitch  angle  of  ship 

DEIN(l) 

t 

Time  of  projectile  in  flight 

DEIN (2) 
DEIN (8) 
DEIN (3) 

X.Y.Z 

Downrange,  cross-range  and  vertical  components 
of  position  of  projectile  in  flight. 

Origin  is  at  the  point  of  firing 

DEIN (4) 
DEIN (9) 
DEIN(5) 

X.Y.Z 

Velocity  components  of  projectile  in  flight 

DEIN(6) 

DEIN(7) 

DC1,DC2 

Drift  integrals 

DEOUT(l) 

• 

t 

Derivative  of  t  with  respect  to  t  (equal  to  1.0) 

DE0UT(2) 

• 

X 

DEIN(4) 

DE0UT(3) 

• 

z 

DEIN(5) 

DEOUT(4) 

DEOUT(9) 

DEOUT(5) 

X.Y.Z 

Components  of  acceleration  of  projectile  in 
flight 

DEOUT(6) 

°C1 

Derivative  of  D^  with  respect  to  t 

DEOUT(7) 

^C2 

Derivative  of  D^  with  respect  to  t 

DE0UT(8) 

Y 

DEIN(9) 
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SYMBOLS  (Continued) 


FORTRAN 

Model 

Text 

DMB 

DMb 

Ship's  speed  perpendicular  to  line  of  fire 

DMBOG 

DMbog 

Cross-range  gun  velocity  component  due  to  gun  throw 

DMHAG 

DMhag 

Gun  velocity  component  perpendicular  to  ship's 
centerline  due  to  gun  throw 

DMHO 

DMho 

Own  ship's  speed 

DMHOG 

DMhog 

Gun  velocity  component  along  ship's  centerline 
due  to  gun  throw 

DMRH 

DMrh 

Own  ship’s  speed  along  line  of  fire 

DMRHOG 

DMrhog 

Horizontal  gun  velocity  component  downrange  due 
to  gun  throw 

DMVOG 

DMvog 

Vertical  gun  velocity  component  due  to  gun  throw 

DRAGC 

DRAGC 

Drag  coefficient  of  projectile 

DRIFT 

DRIFT 

Deflection  perpendicular  to  line  of  fire  due  to 
gyroscopic  effect  on  spinning  projectile 

DS 

DS 

Ballistic  air  density  (percent  of  standard  x  10-^) 

DU 

DU 

Correction  of  projectile  initial  velocity  for  non¬ 
standard  projectile  weight 

DW 

DW 

Variation  of  projectile  weight  from  standard 

DX 

DX 

Horizontal  distance  from  target  future  position 
to  projectile 

DZ 

DZ 

Vertical  distance  from  target  future  position 
to  projectile 

DZO 

VZ(, 

Rate  of  change  of  roll  angle  of  ship 

E 

E 

Stabilized  elevation  angle  of  target  present 
position 

ED 

Ed 

Angle  between  deck  plane  and  line  of  sight. 

measured  in  the  vertical  plane  through  line  of 
sight 
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FORTRAN 

Model 

Text 

SYMBOLS  (Continued) 

EDGP 

Edg' 

Angle  between  deck  plane  and  line  of  fire, 
measured  In  normal  plane  through  line  of  fire 

EG 

Eg 

Gun  elevation  angle  In  stabilized  coordinate  system 

El 

El 

Angle  between  the  horizontal  plane  and  the  deck 
plane,  measured  In  vertical  plane  through  line 
of  sight 

EIO 

Elo 

Pitch  angle  of  ship 

E2 

E2 

Angle  between  horizontal  plane  and  line  to  future 
target  position,  measured  In  vertical  plane  thiough 
line  to  future  target  position 

E4 

E4 

Aiming  position  angle  in  stabilized  coordinate 
system 

G 

G 

Acceleration  due  to  gravity  at  sea  level 

GALT 

GALT 

Geopotential  altitude 

GAMMA 

GAMMA 

Reciprocal  ballistic  coefficient  of  projectile 

IATM 

- 

Atmosphere  type 

ICOR 

- 

Coriolis  type 

LDP 

Ld' 

Angle  between  vertical  plane  through  line  of 
sight  and  normal  plane  through  line  of  fire, 
measured  in  deck  plane  from  vertical  plane  through 
line  of  sight 

LH 

Lh 

Angle  between  line  of  sight  and  line  of  fire 
measured  in  horizontal  plane 

MBG 

Mbg 

Total  linear  deflection  perpendicular  to  line 
of  fire 

MO 

Mode 

Mode  of  operation  of  fire  control  system 

MTHX 

MThx 

Horizontal  displacement  of  target  in  east-west 
direction  during  time  of  flight  of  projectile 
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FORTRAN 

Model 

Text 

MTHY 

MThy 

Horizontal  displacement  of  target  in  north-south 
direction  during  time  of  flight  of  projectile 

MTZ 

MTz 

Vertical  displacement  of  target  during  time  of 
flight  of  projectile 

MXTG 

MxTg 

Horizontal  displacement  of  target  with  respect 
to  gun  in  east-west  direction  during  dead  time 

MYTG 

MyTg 

Horizontal  displacement  of  target  with  respect 
to  gun  in  north-south  direction  during  dead  time 

MZTG 

MzTg 

Vertical  displacement  of  target  with  respect 
to  gun  during  dead  time 

PDA 

Pda 

Displacement  of  gun  from  centerline  of  ship, 
measured  in  deck  plane  perpendicular  to  centerline 

PDO 

Pdo 

Displacement  of  gun  from  director  along  center- 
line  of  ship 

PHA 

Pha 

Component  of  gun  parallax  in  the  horizontal 
plane  perpendicular  to  vertical  plane  through 
own  ship’s  centerline 

PHO 

Pho 

Component  of  gun  parallax  in  the  horizontal  plane 
and  in  the  vertical  plane  through  own  ship's 
centerline 

PTO 

PTO 

Equivalent  density  correction  to  account  for 
variation  of  air  drag  due  to  projectile-fuze 
combinations 

PV 

PV 

Vertical  displacement  of  gun  above  or  below 
director 

PVDP 

Pvd' 

Normal  component  of  gun  displacement  above  deck 

PX 

PX 

East-west  displacement  of  gun  from  director 
in  horizontal  plane 

PY 

PY 

North-south  displacement  of  gun  from  director 
in  horizontal  plane 

QEG 

qEg 

Elevation  spot  correction 

QLH 

qLh 

Deflection  spot  correction 

viii 


FORTRAN 

Model 


Text 


SYMBOLS  (Continued) 


QRHX 

qRhx 

Offset  of  landmark  from  target  in  east-west 
direction 

QRHY 

qRhy 

Offset  of  landmark  from  target  in  north-south 
direction 

QRH4 

qRh4 

Range  spot  correction 

QSTARH 

QSTARH 

Horizontal  distance  of  illuminating  projectile 
beyond  target 

QSTARZ 

QSTARZ 

Vertical  distance  of  illuminating  projectile 
above  target 

QHXS 

Qhxs 

East-west  component  of  QSTARH 

QHYS 

Qhys 

North-south  component  of  QSTARH 

R 

R 

Distance  from  own  ship  to  target  measured 
along  line  of  sight 

RE 

RE 

Radius  of  earth  in  feet 

RH 

Rh 

Projection  of  present  range  in  horizontal 
plane  by  a  vertical  plane  through  line  of  sight 

RHO 

RHO 

Density  of  atmosphere  at  altitude  ZS  or  GALT 

RH2 

Rh2 

Projection  of  R2  in  horizontal  plane  by  a  vertical 
plane  through  R2 

RHX2 

Rhx2 

East-west  component  of  Rh2 

RHY2 

Rhy2 

Nor^h-south  component  of  Rh2 

RH4 

Rh4 

Horizontal  range  to  aim  point 

RV 

RV 

Vertical  range  to  aim  point 

RV2 

RV2 

Height  of  future  position  of  target  above 
horizontal  plane,  measured  in  vertical  plane 
through  R2 

R2 

R2 

Distance  from  gun  to  future  position  of  target 

SYMBOLS  (Continued) 


FORTRAN 

Model 

SR 

T 

TD 

TG 

TK 

TM 

T0L2 

TR 

TS 

T4 

T5 

U 

ucw 

V 

VAA 

W 

WB 

WH 

WHX 


Text 

SR  Slant  range  to  projectile  in  XZ  plane 

T  Air  temperature  (°K)  at  altitude  ZS 

TD  Deviation  of  surface  air  temperature  from 

standard  (°F) 

Tg  Dead  time 

TK  Surface  air  temperature  (°K) 

TM  Surface  air  temperature  (°F) 

T0L2  Upper  limit  of  absolute  value  of  two  successive 

calculated  values  of  time  of  flight  of  projectile 

TR  Air  temperature  (°R)  at  altitude  GALT 

TS  Surface  air  temperature  (°F)  x  10"^ 

T4  Projectile  time  of  flight  at  beginning  of  dead 

time 

T5  Fuze  time  (sec) 

U  Initial  ejection  velocity  of  projectile  from  gun 

UCW  Initial  velocity  of  projectile  corrected  for 

variation  of  projectile  weight  from  standard 

V  Target  speed 

VAA  Velocity  of  projectile  with  respect  to  air 

W  Standard  weight  of  projectile 

Wb  Horizontal  component  of  true  wind  perpendicular  to 

line  of  fire 

Wh  Horizontal  true  wind  speed 

Whx  East-west  component  of  true  wind  in  horizontal 

plane 
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FORTRAN 


Model 

Text 

WHY 

Why 

North-south  component  of  true  wind  in  horizontal 
plane 

WRH 

Wrh 

True  horizontal  ballistic  wind 

X.Y.Z 

x,y,z 

Stabilized  components  of  present  position  of  target 
with  respect  to  director  in  east-west,  north-south 
and  vertical  directions,  respectively 

X.Y.Z 

x,y,z 

Estimated  stabilized  components  of  present  position 
of  target  with  respect  to  director 

XDDOT 

*  *  » 

x,y,z 

Estimated  stabilized  components  of  target  present 

YDDOT 

ZDDOT 

acceleration  with  respect  to  director 

XDOT 

AAA 

*,y,z 

Estimated  stabilized  components  of  target  present 

YDOT 

ZDOT 

velocity  with  respect  to  director 

XDOTTG 

*Tg 

Estimated  stabilized  components  of  target  velocity 

YDOTTG 

yTg 

?.Tg 

with  respect  to  director  at  end  of  dead  time 

ZDOTTG 

XL 

L 

Latitude  of  ship  at  firing  time 

XTG 

xTg 

Stabilized  position  components  of  target  with 

YTG 

yTg 

respect  to  gun  after  dead  time 

ZTG 

zTg 

XW.YW.ZW 

XW,YW,ZW 

Downrange,  cross-range,  and  vertical  components 
of  projectile  velocity  with  respect  to  air 

ZD 

Zd 

Cross  level  angle 

ZO 

Zo 

Roll  ang’e 

ZS 

ZS 

Geometric  altitude  above  sea  level 

Note:  Symbols  used  in  this  report  are,  as  nearly  as  possible, 

consistent  with  standard  fire  control  symbols  published 
in  NAVORD  OP  1700,  Volume  1.0) 
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1 .  INTRODUCTION 


1.1  Background 

The  existing  Gunfire  Control  System  (GFCS)  MARK  68'"'  uses  a 
MARK  47  analog  computer  to  generate  gun  orders  for  the  5-Inch  54-Caliber 
MARK  42  Gun  Mount.  Two  gunfire  control  modes,  air  and  surface,  are 
included  in  the  computer  operation.  The  computer  uses  target  present 
position  data,  ship  and  target  rates  of  motion,  and  ballistic  data  to 
determine  a  line  of  fire  such  that  the  projectile  will  arrive  at  the 
predicted  target  future  position  at  the  end  of  the  projectile  time  of 
flight.  Gun  orders  are  then  computed  and  transmitted  to  the  gun  mount, 
and  fuze-setting  order  is  determined  for  mechanical  time  fuzes.  Under 
the  MARK  68  Gun  Improvement  Program  personnel  at  the  Naval  Surface 
Weapons  Center,  Dahlgren  Laboratory,  worked  toward  converting  the  MARK 
68  Analog  GFCS  to  a  digital  system.  One  result  of  this  work  was  the 
development  of  a  mathematical  model for  digital  gunfire  control  using 
the  NSWC/DL  effective  drag  functions. 

In  order  to  develop  a  more  general  gunfire  control  program,  the 
effective  drag  model  has  been  modified  by  replacing  the  NSWC/DL  functions 
with  fourth-order  Runge-Kutta  numerical  integration  of  three-dimensional 
particle  trajectory  equations^)  for  computation  of  gun  elevation  angle, 
projectile  deflection  including  drift,  and  projectile  time  of  flight. 
Terms  are  included  in  the  particle  trajectory  equations  to  account  for 
effects  of  Coriolis  acceleration,  wind  and  ship's  motion,  and  variation 
of  acceleration  of  gravity  with  altitude.  Drift  is  obtained  by  inte¬ 
gration  of  drift  integrals.  Air  density  and  air  temperature  are  computed 
using  special  equations  for  the  U.S.  Navy  Standard  Atmosphere  (NAST)  or 
the  1962  U.S.  Standard  Atmosphere  (ICAO)  depending  on  which  atmosphere 
is  required  in  the  integration  model.  Projectile  drag  coefficients  for 
NAST  or  ICAO  atmosphere  are  obtained  during  the  integration  process  by 
table  1.  k-up  of  values  stored  in  the  computer  for  the  projectile  of 
interest . 

The  model  in  its  current  form  was  designed  not  for  use  in  a  gun¬ 
fire  control  system,  but  as  a  tool  for  evaluation  of  gunfire  control 
systems.  It  is  expected  that  modifications  will  have  to  be  made  for 
evaluation  of  some  special  features  of  gunfire  control  systems,  and  for 
certain  projectiles  which  do  not  have  free  flight,  are  not  fired  at  the 
predicted  target  future  position,  or  for  which  the  reciprocal  ballistic 
coefficient  is  a  function  of  gun  elevation  angle. 


To  adapt  the  model  for  use  in  a  gunfire  control  system,  it  would 
be  necessary  to  perform  the  ballistic  calculations  with  a  more  efficient 
integration  scheme  than  fourth-order  Runge-Kutta  and  to  functionalize  pro¬ 
jectile  drag  coefficient  using,  for  example,  the  method  described  by  Duke, 
et  al.<5> 

1.2  Objectives 

The  objectives  of  this  report  are  to  present  the  sequence  of 
equations  used  to  solve  gunfire  control  problems  using  numerical  integra¬ 
tion  of  particle  trajectory  equations  and  to  describe  the  method  of  inter¬ 
facing  the  ballistic  equations  with  the  present  position,  prediction,  and 
gun  order  sections  of  the  mathematical  model. 

The  next  section  contains  a  general  description  of  the  modified 
digital  model  and  the  modes  of  operation.  Subsequent  sections  give  the 
sequence  of  calculations  in  the  different  sections  of  the  model,  a  detailed 
description  of  the  ballistic  calculations  using  numerical  integration,  and 
a  discussion  of  the  accuracy  and  limitations  of  the  model. 

Appendix  A  contains  the  FORTRAN  computer  program.  Appendix  B  a 
flow  diagram,  and  Appendix  C  an  input  guide. 


2.  GENERAL  DESCRIPTION  OF  THE  DIGITAL  GUNFIRE  CONTROL  PROGRAM 

2.1  Description 

The  computations  performed  in  this  model  follow  generally  those 
in  Reference  3  except  that  gun  elevation  angle,  projectile  deflection, 
and  projectile  time  of  flight  are  computed  by  numerical  integration  of 
particle  trajectory  equations.  Stabilized  coordinates  of  target  present 
position  are  computed.  The  prediction  section  uses  these  with  target 
velocity,  target  acceleration,  and  projectile  time  of  flight  from  the 
ballistics  section  to  compute  stabilized  coordinates  of  future  position. 
The  ballistics  section  uses  target  future  position  and  ballistic  parame¬ 
ters  to  determine  gun  elevation  angle,  horizontal  deflection  angle  and 
projectile  time  of  flight.  Time  of  flight  is  fed  back  to  the  prediction 
section,  and  the  computation  is  recycled  until  two  successive  values  of 
time  of  flight  are  within  a  specified  tolerance  limit. 

Fuze  time  (T5)  is  taken  as  projectile  time  of  flight  if  dead 
time  was  used  in  determining  future  position. 

Values  of  gun  elevation  angle  and  horizontal  deflection  angle 
computed  at  the  beginning  of  dead  time  are  input  to  the  gun  order 
section  where  gun  train  order  and  gun  elevation  order  are  computed. 

2.2  Modes  of  Operation 

Since  the  ballistic  calculations  are  accomplished  using  numerical 
integration,  fire  control  problems  can  be  solved  for  each  mode  of  firing 
any  standard  projectile  for  which  drag  data  and  the  reciprocal  ballistic 
coefficient  are  available.  No  distinction  is  made  between  air  and 
surface  modes.  Reduced  charge,  full  charge  and  super  charge  computations 
can  be  made  by  using  corresponding  drag  data  and  proper  values  for 
reciprocal  ballistic  coefficient. 

Program  branches  are  provided  for  both  reverse  slope  and  star 
shell  modes.  It  should  be  noted  that  projectile  drag  characteristics 
may  vary  between  normal  trajectory  and  high  angle  firings.  These 
variations  will  have  to  be  provided  by  the  user.  The  star  shell  geo¬ 
metry  is  described  in  detail;  however,  implementation  of  drag  variation 
after  formation  of  the  star  must  be  provided  also  by  the  user. 

Some  modifications  would  be  required  for  firing  a  guided  pro¬ 
jectile,  the  flight  of  which  is  affected  by  thrust  and  guidance. 
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3.  PRESENT  POSITION  SECTION 

In  the  present  position  section,  data  measurements  from  the  deck 
instruments  are  used  to  compute  stabilized  polar  coordinates  of  target 
present  position  with  respect  to  the  director;  stabilized  components  of 
gun  displacement  from  director  along  and  perpendicular  to  the  ship's 
centerline  and  in  the  east-west,  north-south  and  vertical  directions; 
and  Cartesian  coordinates  of  target  present  position  with  respect  to  the 
director. 

Provisions  are  made  for  offset  and  range  spot  corrections  in  the 
surface  mode.  Projectile  initial  velocity  and  ballistic  air  density  are 
modified  to  compensate  for  variation  of  projectile  weight  from  standard. 
Ballistic  air  density  is  modified  to  compensate  for  special  projectile- 
fuze  combinations.  East-west  and  north-south  components  of  true  wind, 
and  roll  and  pitch  angles  of  the  ship  are  computed. 

Cartesian  coordinates  of  present  position  are  computed  with  respect 
to  the  director  at  a  given  time  and  used  as  best  estimates  of  target 
present  position  along  with  velocity  and  acceleration  coordinates  from 
whatever  source  is  available  to  the  user.  These  estimated  values  and 
appropriate  parallax  components  are  used  to  determine  coordinates  of 
target  position  with  respect  to  gun  at  end  of  dead  time. 

Input  to  the  present  position  section  includes  range  to  present 
position  of  target  (R),  director  train  (Bd),  director  elevation  (Ed), 
level  angle  (Ei),  cross-level  angle  (Zd),  heading  of  ship  (Cqo) ,  ship's 
course  (Co),  ship's  speed  (DMho) ,  dead  time  (Tg),  east-west  offset 
(qRhx) ,  north-south  offset  (qRhy) ,  range  spot  correction  (qRh4) ,  gun 
parallax  displacements  (Pdo,  Pda,  and  Pvd'),  projectile  initial  velocity 
(U),  standard  weight  of  projectile  (W),  variation  of  projectile  weight 
from  standard  (DW),  ballistic  air  density  (DS),  equivalent  density 
correction  for  particular  projectile-fuze  combinations  (PTO),  true  wind 
speed  (Wh)  and  direction  (Bwy) ,  and  constant  AM  used  in  correction  of 
projectile  initial  velocity  for  variation  of  projectile  weight  from 
standard.  Horizontal  distances  of  illuminating  projectile  from  target 
in  the  east-west  (Qhxs)  and  north-south  (Qhys)  directions,  and  distance 
of  illuminating  projectile  above  target  (QSTARZ)  are  set  to  zero  for 
applications  other  than  star  shell. 

The  sequence  of  computations  performed  in  the  present  position 
section  follows.  Symbols  are  defined  in  a  separate  section  of  this 
report  (see  Contents)  and  are  generally  the  same  as  those  in  Reference  3. 


Eio  =  sin-1 (sin  Ei  •  cos  Bd  +  cos  Ei  •  sin  Bd  •  sin  Zd) 

Zo  =  sin_1[(cos  Bd  •  cos  Ei  •  9in  Zd  -  sin  Ei  •  sin  Bd)/cos  Eio] 

B  =  tan-1 
By  =  Cqo  +  B 
E  =  Ed  -  Ei 
Rh  =  R  •  cos  E 
x  =  Rh  •  sin  By  +  qRhx 

y  -  Rh  •  cos  By  +  qRhy 

By  =  tan-1  (x/y) 

B  =  By  -  Cqo 

Rh  =  (x2  +  y2)^  +  qRh4 

Zd  =  sin_1[sin  B  •  cos  Zo  •  sin  Eio  +  cos  B  •  sin  Zo] 

Ei  =  sin_1[(cos  B  •  cos  Zo  •  sin  Eio  -  sin  B  •  sin  Zo)/cos  Zd] 

Bd  =  tan-1[(cos  Ei)/(cot  B  •  cos  Zd  +  sin  Ei  •  sin  Zd)  ] 

Pho  =  Pdo  •  cos  Eio  +  Pvd 1  •  sin  Eio  •  cos  Zo  +  Pda  •  sin  Eio  •  sin  Zo 
Pha  =  Pda  •  cos  Zo  -  Pvd'  •  sin  Zo 
PX  =  Pho  •  sin  Cqo  +  Pha  •  cos  Cqo 

PY  =  Pho  ■  cos  Cqo  -  Pha  •  sin  Cqo 

PV  =  Pvd'  •  cos  Eio  •  cos  Zo  +  Pda  •  sin  Zo  •  cos  Eio  -  Pdo  •  sin  Eio 
x  =  Rh  •  sin  By 
y  =  Rh  '  cos  By 
z  =  R  •  sin  E 

x,  y  and  z  are  stabilized  Cartesian  present  position  coordinates 
01  the  target  with  respect  to  the  director  at  a  given  time.  These  are  ingu£  as 
best  estimates  of  fca^get  present  position  (x,y,z)  along  with  velocity  (x  y  z) 
and  acceleration  (x,y,z)  coordinates  from  whatever  source  is  available  to  the 
user . 
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•  \rx2  ~rj 


DU  *  -AM  •  DW  •  U/W 


DDS  =  -DS  •  DW/W 


UCW  =  U  +  DU 


DS  =  DS  +  DDS  +  PTO/ 100 
Whx  =  -Wh  ■  sin  Bwy 

Why  ■  -Wh  •  cos  Bwy 

b  ,  A  _ 

MxTg  =  x  •  Tg  +  $  •  x  *  Tg 

MyTg  =  y  •  Tg  +  |  .  y  *  Tg2 

A  A 

MzTg  =  z  •  Tg  +  £  •  z  •  Tg3 

A 

xTg  =  x  +  MxTg  -  PX 

yTg  *  y  +  JfyTg  -  PY 
A 

zTg  =  z  +  MzTg  -  PV 
ByTg  =  tan-1  (xTg/yTg) 


4.  PREDICTION  SECTION 

The  future  position  of  the  target  is  the  target  position  with  respect 
to  the  gun  at  the  end  of  the  time  of  flight  of  the  projectile.  It  is  ob¬ 
tained  by  adding  the  displacement  of  the  target  during  the  time  of  flight 
to  the  position  of  the  target  with  respect  to  the  gun  at  the  end  of  dead 
time.  In  normal  modes  of  operation  (i.e.,  other  than  star  shell),  the  pro¬ 
jectile  is  delivered  to  the  future  position  of  the  target.  In  the  star 
shell  mode,  the  projectile  is  delivered  to  some  point  beyond  and  above 
the  threat  depending  upon  the  projectile  being  used  and  the  desired  result. 

The  sequence  of  computations  performed  in  the  prediction  section 
follows . 


=  (xS  +  ya 


+  I)* 


T4  =  R/(U  +  V) 

A  A  A 

xTg  =  x  +  x  •  Tg 

A  A  A 

yTg  =  y  +  y  •  Tg 

A  A  A 

•  m  •  •• 

zTg  =  z  +  z  •  Tg 

A  A 

MThx  =  iTg  •  T4  +  •  x  •  T42  +  DMho  •  T4  •  sin  Co 


MThy  =  yTg  •  T4  +  £ 

A 

MTz  =  zTg  •  T4  +  i 


A 

y 


T42  +  DMho  •  T4  *  cos  Co 
T43 


star 

shell,  recompute  MThx  and 

MThy, 

and  compute 

Qhxs 

and  Qhys 

MThx 

A 

=  (xTg  +  DMho  •  sin 

Co)  * 

(T4  +30) 

+  $•£•• 

(T4 

+  30)s  - 

MThy 

A 

=  (yTg  +  DMho  •  cos 

Co)  • 

(T4  +30) 

+  ?  •  y  • 

(T4 

+  30) 2  - 

x' 

=  xTg  +  (30  +  T4)  • 

xTg  + 

A 

i  •  x 

• 

(30  +  T4)2 

y ' 

=  yTg  +  (30  +  T4)  • 

A 

yTg  + 

A 

£  ’  y 

• 

(30  +  T4)2 

R' 

=  (x,a  +  y'2)^ 

Whx  •  30 


QSTARZ  =  2000,  QSTARH  =  5400. 

(Values  of  QSTARZ  and  QSTARH  may  be  varied  as  desired) 
Qhxs  =  QSTARH  •  x'/R' 

Qhys  =  QSTARH  •  y'/R' 
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The  star  shell  equations  above  are  for  the  5-lnch  54  Caliber 
MARK  48  Illuminating  Projectile;  introduction  of  other  rounds  may 
require  a  different  sequence  of  equations. 

The  equations  below  are  for  all  projectiles  and  modes  of  firing. 
For  projectiles  other  than  star  shell,  Qhxs  ■  Qhys  •  QSTARZ  -  0. 


Rhx2  *  xTg  +  MThx  +  Qhxs 
Rhy2  -  yTg  +  MThy  +  Qhys 
RV2  =  zTg  +  MTz  +  QSTARZ 
By2  *  tan-* (Rhx2/Rhy2) 

R2  =  [(Rhx2)2  +  (Rhy2) 2  +  (RV2)2]*5 
E2  =  sin-1 (RV2/R2) 


5 .  BALLISTICS  SECTION 

The  ballistics  section  uses  future  position  coordinates,  projectile 
constants,  projectile  drag  coefficients  and  other  input  quantities  shown 
below  to  calculate  gun  bearing,  horizontal  deflection  angle,  gun  eleva¬ 
tion  angle  and  projectile  time  of  flight.  Gun  elevation  angle,  projec¬ 
tile  time  of  flight,  and  projectile  deflection  are  calculated  by  fourth- 
order  Runge-Kutta  integration  of  three-dimensional  particle  trajectory 
equations.  Time  of  flight  is  fed  back  to  the  prediction  section  and 
calculations  are  repeated  until  two  successive  values  of  time  of  flight 
are  within  a  specified  tolerance.  When  this  condition  is  satisfied, 
time  of  flight  is  taken  as  fuze  time.  The  calculation  is  repeated 
for  the  present  position  of  the  ship,  using  Tg  equal  to  zero,  to  obtain 
values  for  calculation  of  gun  orders. 

Output  from  the  ballistics  section  includes  horizontal  deflection 
angle  (Lh)  and  gun  elevation  angle  (Eg)  to  the  gun  order  section,  and 
time  of  flight  (T4)  to  the  prediction  section. 

5.1  Input  Not  Available  From  Present  Position  and  Prediction  Sections 

Input  to  the  ballistic  section  that  is  not  generated  in  the  present 
position  and  prediction  sections  includes  air  temperature  (TS)  ,  earth 
constants  (L,G),  tolerance  used  in  obtaining  solutions  of  gun  train 
angle  and  time  of  flight  of  projectile  (T0L2),  projectile  drift  con¬ 
stant  (CANS),  and  rates  of  pitching  (DEio) ,  rolling  (DZo) ,  and  yawing 
(DCqo)  of  the  ship. 

RV  is  set  equal  to  RV2  and  the  initial  value  of  Bgy  is  taken  as 

By2 . 


Horizontal  range  to  aim  point  and  aiming  position  angle  are 

given  by 


Rh4  =  Rh2  •  cos (Bgy  -  By2) 


and 

E4  =  tan_1(RV/RhA), 


respectively. 

5.2  Gun  Velocity  Components  Due  to  Gun  Throw 

Gun  throw  is  the  linear  velocity  of  the  gun  resulting  from 
rolling,  pitching  and  yawing  of  the  ship.  The  components  of  this 
velocity  along  and  perpendicular  to  the  centerline  of  the  ship 
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are 


and 


respectively , 


DMhog  =  PV  •  DEio  -  Pha  •  DCqo 


DMhag  =  -PV  •  DZo  +  Pho  •  DCqo, 


The  components  downrange,  cross-range  and  vertical,  due  to  gun 
throw,  are 

DMrhog=  DMhog  •  cos  Bg  +  DMhag  •  sin  Bg, 

DMbog  =  DMhag  ■  cos  Bg  -  DMhog  •  sin  Bg, 
and 

DMvog  =  -Pho  •  DEio  +  Pha  ’  DZo, 
respectively,  where 

Bg  =  Bgy  -  Cqo . 

Initial  velocity  is  recalculated  to  correct  for  the  effects  of  gun  throw 
using 

U  =  [(UCW  •  cos  Eg  +  DMrhog) 2  +  (UCW  •  sin  Eg  +  DMvog)*]^. 
5.3  Wind  and  Ship  Velocity  Components 

Components  of  true  wind,  downrange  and  cross-range,  are 
Wrh  =  Why  *  cos  Bgy  +  Whx  •  sin  Bgy 

and 

Wb  =  —Why  •  sin  Bgy  +  Whx  •  cos  Bgy, 

respectively . 

Components  of  the  velocity  of  the  ship  downrange  and  cross-range 


are 


and 


DMrh  =  DMho  ‘  cos  (Co  -  Bgy) 


DMb  =  DMho  ’  sin  (Co  -  Bgy) , 


respective ly . 
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5.4  Gun  Elevation  Angle,  Projectile  Deflection,  and  Projectile 
Time  of  Flight 

The  slant  range  to  the  aim  point  is 
SR4  =  (Rh42  +  RV2)*2 

An  approximate  value  of  gun  elevation  angle  is  given  by 
Eg  =  tan_1(RV/Rh4)  +0.5  sin-1 (ANG) 


where 


ANG  =  64.4  ■  Rh4/U2  . 

If  ANG  >  1.,  set  ANG  =  1.0.  For  high  angle  or  reverse  slope  firing. 

Eg  =  7r/2  -  Eg. 

Initialize  t  =  1.,  t=0.,  X  =  0 . ,  Y  =  0 . ,  Z=0.,  £  =  U  •  cos  Eg  + 
DMrh,  Y  =  DMb  +  DMbog,  t  =  U  •  sin  Eg,  Dcl  =  0.,  and  Dc2  =  0.  X,  Y,  Z  and 
£,  tf,  i  are  position  and  velocity  components  of  the  projectile,  respectively 
X  is  positive  down  range,  Y  is  positive  to  the  right  as  seen  from  the  point 
of  firing,  and  Z  is  positive  upward. 

The  downrange,  cross-range,  and  vertical  components  of  velocity  of 
the  projectile  relative  to  the  air  are 

XW  =  i  -  Wrh,  YW  =  Y  -  Wb 

and  ZW  =  Z,  respectively.  The  velocity  of  the  projectile  relative  to  the 
air  is 

VA  =  [(XW)8  +  (YW)8  +  (ZW)2]^. 

Air  temperature,  air  density,  and  velocity  of  sound  are  required  at 
each  step  in  the  numerical  integration  of  the  trajectory  equations.  These 
are  calculated  for  the  desired  atmosphere  (NAST  or  ICAO)  using  equations 
given  below. 

NAST  Atmosphere 


Surface  air  temperature  is  given  in  degrees  Kelvin  by 


TK  =  5  •  TM/9  +  255.23 


where  TM  =  100  ■  TS  . 
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Temperature  and  velocity  of  sound  at  geometric  altitude  ZS  above 
sea  level  are  given  by 

T  =  TK  -  0.001982  •  ZS 
and 

CS  =  CSS  *  (T/288)5* , 
respectively,  where 

ZS  =  Z  +  X3/ (41812000) . 

Z  is  geometric  altitude  above  the  tangent  plane  at  the  point  of  firing 
the  projectile. 

Density  of  the  atmosphere  at  altitude  ZS  is  given  in  terms  of 
NAST  ballistic  air  density  (DS)  by 

RHO  =  0.07513  *  EXP  (-0.0000031582-ZS)  •  DS . 

ICAO  Atmosphere 

Formulas  for  temperature,  pressure,  and  density  of  the  ICAO 
atmosphere  are  given  in  Reference  6.  Constants  were  determined  by  Loren 
J.  McAnelly  (NSWC/DL)  for  a  modified  version  of  these,  and  are  given  below 
for  several  ranges  of  altitudes.  Symbols  have  been  changed  to  conform  to 
usage  in  this  report . 

Values  are  given  in  terms  of  geopotential  altitude  (GALT)  where 
GALT  =  RE  •  ZS/  (RE  +  ZS)  .  RE  is  radius  of  the  earth,  20855531  feet, 
and  ZS  =  Z  +  X2 / (41812000) . 

Standard  temperature  (TR)  and  standard  pressure  (P)  at  altitude 
GALT  are  computed  as  follows  for  different  ranges  of  geopotential  alti¬ 
tude  : 


-5000  £  GALT  <  36089  feet 

TR  =  518.67  -  0.00356616  •  GALT 
P  =  1013.25  (518 .67/TR) "B  -a56877 
36089  5  GALT  <  65616.8  feet 
TR  =  389.97 

P  =  226.32  EXP  [-0.0000480634 (GALT -36089. 239)] 
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65616.8  ^  CALT'c  104986  feet 

TR  =  389.97  +  0.00054864  (GALT-65616.789) 

P  =  54.7487  (389.97/TR)34- 1631  9 
TR  is  in  degrees  Rankine  and  P  is  in  millibars. 

Standard  air  density  for  the  geopotential  altitude  range 

-5000  £  GALT  <  104986  feet  is  given,  in  pounds  per  cubic  foot, 
by 

RHO  =  0.0391462  (P/TR) . 

For  nonstandard  ICAO  ballistic  density  (DS) 

RHO  =  RHO  •  DS. 

Velocity  of  sound  at  altitude  GALT  is  given  by 
CS  =  1116.49  [(TR  +  TD)/ (5 18. 67)]  % 

where  TD  is  deviation  of  surface  air  temperature  from  standard  in  degrees 
Fahrenheit.  Mach  number  is  computed  for  either  NAST  or  ICAO  atmosphere 
using 

CM  =  VA/CS. 

The  projectile  drag  coefficient  (DRAGC)  is  obtained  as  a  function 
of  Mach  number,  for  the  atmosphere  being  used,  by  table  look-up  of 
values  stored  in  the  computer. 

The  equations  of  motion  of  the  projectile  in  three  dimensions  are 
X  =  CONS  •  XW  -  (AY  •  Z  -  AZ  •  Y)  -G  *X/RE, 

Y  =  CONS  •  YW  -  (AZ  •  X  -  AX  •  Z)  -  G  •  Y/RE, 

and 

Z  =  CONS  •  ZW  -  (AX  •  Y  -  AY  -X)  -  G(1  -  2Z/RE) 

where 

AX  =  -0.00014584  •  cos  L  •  cos  Bgy, 

AY  =  +0.00014584  •  cos  L  •  sin  Bgy, 

and 

AZ  =  -0.00014584  •  sin  L. 

RE  =  20855531  feet  and 
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CONS  =  -  tt/8  •  GAMMA  •  RHO  •  VA  *  DRAGC/144 . 

GAtttlA  is  the  reciprocal  ballistic  coefficient  of  the  projectile. 

Projectile  drift  is  obtained  from  the  solutions  of  the  equations 
Doi  =  l/VA3  and  D68  =  X/VA2. 

Integration  of  the  trajectory  and  drift  equations  to  X  =  Rh4  (Z  =  RV 
for  reverse  slope)  yields  values  of  X,  Y,  i,  Del,  Doa,  X,  Y,  Z  and  t.  The 
integration  method  used  is  standard  fourth-order  Runge-Kutta  with  a  fixed 
integration  interval. 

Drift  is  computed  using 

DRIFT  =  CANS  (X  *  Dal  -  Daa) 

where  CANS  is  the  drift  constant  for  the  projectile . 

Let  DZ  =  RV  -  Z  (DX  =  Rh4  -  X  for  reverse  slope)  and  SR  =  (X2  +  Z2)^ 

then 

Eg  =  Eg  +  tan_1[DZ  *  (cos  E4)/SR], 

(For  reverse  slope  Eg  =  Eg  +  tan_1[DX  *  (cos  Eg)/Rh4].) 

The  integration  is  repeated  using  the  new  value  of  Eg  until  DZ 
(DX  for  reverse  slope)  is  less  than  some  preassigned  value  (say  0.1  ft). 

When  this  condition  is  satisfied,  set  T4  =  t  and  compute  deflection  of  the 
projectile  perpendicular  to  the  line  of  fire  using 

Mbg  =  Y  +  DRIFT . 
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5.5  Computation  of  Train  Angle  and  Horizontal  Deflection  Angle 

The  gun  train  angle  in  the  horizontal  plane  is  the  azimuth  of 
the  line  of  fire  measured  clockwise  from  north.  It  is  the  sum  of  the 
azimuth  of  the  target  future  position  (By2)  and  the  negative  of  the 
angle  (DBy)  required  to  compensate  for  deflection  of  the  projectile  in 
flight.  The  train  angle  is  given  by 

Bgy  =  -DBy  +  By2 


where 


DBy  =  tan-1 (Mbg/Rh4) 

When  the  train  angle  is  computed,  the  horizontal  deflection 
angle  (Lh)  is  obtained  from  the  equation 

Lh  =  Bgy  -  ByTg 

Lh  is  used  in  computation  of  gun  orders. 

Time  of  flight  is  fed  back  to  the  prediction  section  and  the 
computation  loop  is  recycled  until  two  successive  values  of  time  of 
flight  are  within  a  given  tolerance  (T0L2).  When  this  condition  is  met, 
T5  is  set  equal  to  T4  if  future  position  was  predicted  at  the  end  of 
dead  time.  If  future  position  was  predicted  at  the  beginning  of  dead 
time  using  Tg  =  0,  Eg  and  Lh  are  output  to  the  gun  order  section  and  gun 
orders  are  computed. 

The  method  of  closure  of  the  solution  for  time  of  flight  of  the 
projectile  is  descrioed  in  Reference  3,  Appendix  J. 
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6.  GUN  ORDER  S£CTION 

6.1  Fuze  Time 

When  future  position  is  predicted  at  the  end  of  dead  time,  the 
projectile  flight  time  (T4)  is  taken  as  fuze  time  (T5)  . 

6.2  Gun  Orders 

Lh  and  Eg  are  corrected  by  adding  deflection  (qLh)  and  elevation 
(qEg)  spot  corrections,  if  available. 

Then 


and 


Lh  =  Lh  +  qLh 


Eg  =  Eg  +  qEg. 

Gun  elevation  order  is  computed  by  substituting  values  of  Lh  and 
Eg  into  the  equation 

Edg'  =  sin-1  [-sin  Zd  •  sin  Lh  *  cos  Eg  +  cos  Zd  • 

(sin  Ei  •  cos  Eg  •  cos  Lh  +  sin  Eg  •  cios  Ei)]. 


To  obtain  gun  train  order,  compute 


Ld' 


_1 [cos  Zd-sin  Lh‘Cos  Eg+sin  Zd*(sin  Ei»cos  Eg’cosLh+sinEg-cosEi)  ; 
tan  !.  cos  Eg  *  cos  Ei  •  cos  Lh  -  sin  Eg  •  sin  Ei  J 


then 


Bdg 1  =  Ld'  +  Bd. 
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7.  ACCURACY  OF  MODEL 


Accuracy  of  the  numerical  integration  gunfire  control  model  has  been 
checked  for  static  targets  against  range  tables  and  standard  particle 
models,  and  for  moving  targets  against  the  effective  drag  model^)  and 
OP  3729(2'  results. 

Typical  results  of  the  range  table  comparison  are  shown  below  for 
standard  conditions. 

Time  of  Flight  (Sec)  Gun  Elevation  Angle  (Min) 


Numerical 

Numerical 

Range 

Range 

Integration 

Range 

Integration 

(Ft) 

Table 

Model 

Table 

Model 

6000 

2.55 

2.549 

57.5 

57.6 

12000 

5.44 

5.440 

125.6 

125.6 

18000 

8.75 

8.749 

207.2 

207.3 

24000 

12.57 

12.571 

306.4 

306.6 

36000 

22.23 

22.235 

581.0 

581.3 

42000 

28.15 

28.154 

770.2 

770.6 

45000 

31.30 

31.304 

879.1 

879.5 

Comparison  of  the  model  with  a  standard  particle  model 

is  shown  below  for 

nonstandard  conditions. 

Time  of  Flight  (Sec) 

Gun  Elevation  Angle  (Min) 

Position 

Numerical 

Numerical 

Range 

Angle 

Particle  Integration 

Particle 

Integration 

(Yds) 

(Deg) 

Model  Model 

Model 

Model 

5527.0 

78.87 

8.11  8.114 

4768.8 

4769.2 

6482.7 

34.23 

9.90  8.898 

2253.6 

2253.8 

7678.7 

7.52 

11.10  11.099 

704.4 

704.0 

11890.7 

2.61 

19.39  19.393 

619.2 

619. 1 

These  tabulations  demonstrate  that  the  numerical  integration  gunfire 
control  model  duplicates  results  given  in  the  range  table  and  those  com¬ 
puted  using  a  standard  particle  model. 

In  NAVORD  OP  3729^  problems  are  provided  for  A-Tests  for  the 
Gunfire  Control  System  MARK  68  MOD  13  firing  the  5-Inch  projectile 
MARK  41  with  VT  fuze  MARK  73.  Among  these  are  twelve  AA-800  (target 
speeds  to  800  knots)  and  eight  AA-2000  (target  speeds  to  2000  knots) 
problems  numbered  11-22  and  26  -  33,  respectively.  For  each  problem 
own  ship  data,  target  data,  atmospheric  data  and  ballistic  corrections 
are  given.  Solutions  are  given  for  each  problem  for  fuze  setting  (T5)  , 
gun  train  order  (BDG').  and  gun  elevation  order  (Edg1). 


17 


The  cwenty  AA  problems  were  solved,  without  Coriolis  effects,  by 
the  mathematical  model  described  in  this  report.  Solutions  were  obtained 
for  no  Coriolis  acceleration  because  solutions  given  in  OP  3729  neglected 
Coriolis  effects.  For  each  problem  T5,  Bdg'  and  Edg'  were  computed. 

Table  1  compares  these  solutions  with  those  given  in  OP  3729  and  with 
those  obtained  using  the  mathematical  model  described  in  Reference  3. 

The  columns  of  the  table  show  for  each  problem  the  deviation  of  the  OP 
results  and  the  results  obtained  with  the  digital  model (3)  from  the 
results  obtained  with  this  model  which  solves  ballistics  by  numerical 
integration  of  particle  trajectory  equations.  Root-mean-square  and 
average  values  of  the  differences  are  given  at  the  bottom  of  each 
column. 

The  root-mean-square  values  for  deviation  of  the  digital  model  from 
results  of  numerical  integration  are  about  twenty  to  forty  percent  of 
those  for  deviation  of  the  OP  results  from  numerical  integration.  This 
Is  consistent  with  conclusions  in  Reference  3.  Average  values  of  the 
deviations  for  the  digital  model  are  less  than  half  those  for  the  OP. 

Results  in  Table  1  indicate  close  agreement  between  the  digital 
model^)  and  the  numerical  integration  model.  In  view  of  this  agreement, 
it  can  be  safely  concluded  that  the  numerical  integration  model  computes 
the  correct  gun  orders  for  the  A-Test  problems. 

Additional  checks  on  the  accuracy  of  the  model  including  computation 
of  Coriolis  effects  are  reported  in  Reference  7. 


18 


TABLE  1 


DEVIATIONS  OF  OP  AND  DIGITAL  MODEL  SOLUTIONS  OF 
A-TEST,  AA  MODE  GUNFIRE  CONTROL  PROBLEMS  FROM 
SOLUTIONS  OBTAINED  BY  NUMERICAL  INTEGRATION 


(No  Coriolis) 


T5(Sec) 

BDG'(Min) 

EDG'(Min) 

A-Test 

Deviation 

From 

Deviation 

From 

Deviation  From 

Problem 

Numerical  Integration 
Digital 

Numerical  Integration 
Digital 

Numerical  Integration 
Digital 

OP  3729 

Model 

OP  3729 

Model 

OP  3729 

Model 

1 1 

-.001 

.004 

-.1 

.3 

1.5 

-> 

12 

.028 

-.001 

3.3 

.5 

-3.5 

-.5 

13 

.039 

-.003 

1.7 

.3 

.8 

.0 

14 

.020 

.000 

.0 

.2 

-.1 

15 

-.002 

-.008 

-1.2 

2.9 

-.7 

-.9 

16 

.008 

-.006 

-.6 

.4 

1.4 

-.7 

17 

.015 

-.008 

.0 

.4 

1.4 

-.5 

18 

.005 

-.006 

.3 

-.6 

.0 

-.1 

19 

-.008 

.006 

.1 

.0 

-.2 

-.1 

20 

-.001 

.008 

.1 

-.5 

.3 

-.1 

21 

.021 

-.002 

-1.8 

.8 

.8 

.0 

22 

.013 

-.009 

-1.2 

2 

1.2 

-.7 

26 

.029 

-.002 

.5 

-.9 

.1 

.1 

27 

.028 

-.003 

3.8 

-1.1 

-.4 

-1.1 

28 

-.001 

.013 

-1.8 

-2.0 

.3 

-.5 

29 

.004 

-.003 

-1.3 

-.4 

1.1 

-.1 

30 

.043 

.007 

-.7 

-1.1 

.9 

-.4 

31 

.028 

.020 

-7.9 

-.4 

7.2 

3.7 

32 

-.009 

-.002 

-5.9 

-.1 

-  2  2 

_  i 

33 

.021 

-.005 

.5 

“> 

1.0 

T 

RMS 

.021 

.007 

2.6 

1.0 

3.9 

.9 

AV* 

.016 

.006 

1.6 

.7 

1.3 

.5 

Average  of  absolute  values. 
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8.  DISCUSSION 


With  this  mathematical  model,  problems  may  be  solved  using  either 
NAST  or  ICAO  atmosphere.  It  is  necessary  to  use  drag  data  and  ballistic 
density  corresponding  to  the  chosen  atmosphere. 

To  compute  gun  orders  for  star  shell  and  guided  projectiles,  the 
numerical  integration  model  would  have  to  be  modified.  Modifications  of 
the  model  would  be  required,  also,  to  check  special  characteristics  of 
some  gunfire  control  systems. 

For  some  projectiles,  the  reciprocal  ballistic  coefficient  is  a 
function  of  gun  elevation  angle.  A  program  change  would  be  required  to 
compute  gun  orders  for  such  a  projectile. 

Prior  to  the  writing  of  this  report,  the  numerical  integration  gun¬ 
fire  control  model  described  herein  was  used,  with  excellent  results,  to 
predict  gun  orders  for  evaluation  of  the  Coastal  Patrol  and  Inter¬ 
diction  Craft  (CPIC)  gunfire  control  system.  It  is  felt  that  this  working 
validation  and  the  results  reported  in  Section  7  constitute  an 
acceptable  checkout  for  the  model  to  allow  it  to  be  used  as  a  system 
design  tool. 
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INPUT  BALLISTIC  QUANTITIES, 
CONTROL  AND  COMPARISON  DATA 


FUTURE  POSITION  CALCULATION 


BALLISTICS 


CALCULATE  6UN  ORDERS 


INPUT  GUIDE 


Card 


Group 

Card 

Column 

Format 

Symbol 

Description 

1 

1 

1-10 

A4 

PNAME 

Projectile  type 

2 

1-6 

F6.4 

CANS 

Projectile  drift  constant 

7-12 

F6.4 

GAMMA 

Reciprocal  ballistic  coefficient 

13-18 

F6.4 

G 

Sea  level  acceleration  of  gravity 
(ft/secJ ) 

19-24 

F6.4 

CSS 

Sea  level  velocity  of  sound 
(ft/sec) 

25-30 

F6.4 

AM 

Projectile  weight  correction 
factor 

31-36 

F6.4 

PTO 

Projectile-fuze  correction  factor 
(%) 

2 

1 

1-5 

15 

NX 

Number  of  ig  coefficient  pairs 

3 

1 

1-72 

12F6.4 

VAA(l) 

Mach  numbers  in  ascending 

2 

a 

• 

a 

order  (NX  values,  12  per  card) 

4 

a 

1 

a 

1-72 

a 

12F6.4 

DRAG(I) 

Drag  coefficient  (CD) 
corresponding  to  VAA(l)  (NX 
values,  1 2  per  card) 

5 

1 

1-5 

15 

MO 

Mode:  1,  AA;  3,  surface, 
regular  charge;  4,  surface, 
reduced  charge;  10,  MK  48  WP 

2 

1-5 

15 

I  ATM 

Type  of  atmosphere:  1 .  NAST; 
2,  ICAO 

3 

1-5 

15 

ICOR 

Coriolis  type:  1,  none:  2,  with 

4 

1-5 

15 

IPRNT 

Print  type:  0,  most;  1 ,  more; 2, 
least;  5.  table 

5 

1-5 

15 

1DID 

Diddle  type:  1.  without;  2. 
with*;  3,  with 

6 

1-5 

15 

NUMB 

Number  of  problems  to  be  run 

6 

1 

1-5 

15 

IDEN 

Identification  number  of 
problem 

Not  used  in  any  calculations.  Used  only  in  print  out. 


C-l 


Card 


Card 

Column 

Format 

Symbol 

Description 

2 

1-10 

F10.3 

BO 

Director  bearing  (deg) 

11-20 

FI  0.3 

BWY 

Wind  direction  (deg) 

21-30 

F10.3 

CO 

Ship  course  (deg) 

3140 

F10.3 

DMHO 

Ship  speed  (kt) 

41-50 

F10.3 

Ed 

Director  elevation  (deg) 

51-60 

FI  0.3 

Ei 

Level  angle  (deg) 

61-70 

F10.3 

IV 

Initial  velocity  (ft/sec) 

3 

1-10 

F10.3 

PDO 

Gun  parallax  along  ship’s 
centerline  (ft% 

11-20 

F10J3 

PT 

Projectile  typo  (%) 

21-30 

FI  0.3 

PVDP 

Vertical  parallax  (ft) 

3140 

F10.3 

R 

Target  range  (yd) 

41-50 

F10.3 

RHO 

Air  density  (%) 

51-60 

F10.3 

TEM 

Air  temperature  (°F) 

61-70 

F10.3 

TG 

Dead  time  (sec) 

4 

1-10 

F10.3 

W 

Projectile  weight  (lb) 

11-20 

F10.3 

WH 

Wind  speed  (kt) 

21-30 

F10.3 

ZD 

Cross  level  (deg) 

3140 

F10.3 

DMV 

Vertical  target  velocity  (ft/sec) 

41-50 

F10.3 

DMhx 

East-west  target  velocity  (ft/sec) 

51-60 

F10.3 

DMhy 

North-south  target  velocity 
(ft/sec) 

61-70 

F10.3 

T50P 

Fuze  time  (sec)** 

5 

1-10 

F10.3 

EDGPOP 

Elevation  gun  order  (min)** 

11-20 

F10.3 

BDGPOP 

Bearing  gun  order  (min)** 

The  user  now  repeats  the  5  cards  of  card  group  6  until  “NUMB”  input  sets  have  been 

provided. 


“These  three  input  values  will  be  compared  to  the  actual  results  of  the  gunfire  control  solution  in  the  IPRINT  =  2  print 
option.  They  do  not  alter  the  calculations. 


C-2 


